Absence of topological degeneracy in the Hubbard model on honeycomb lattice 



O 

(N 

o 
O 



I 

o 
o 



> 

(N 
(N 
(N 

\6 
o 
o 



Tao Li 

Department of Physics, Renmin University of China, Beijing 100872, P.R.China 

(Dated: October 26, 2010) 

It is shown that the unique sign structure of the ground state of the Hubbard model on honeycomb 
lattice, which is shown to be insensitive to the trapped Z2 gauge flux when the system is defined 
on a torus, may cause the absence of topological degeneracy on this bipartite system. Examples of 
variational Mott insulating state on the honeycomb lattice are given to illustrate the close relation 
between the sign structure of the ground state and the (absence of) topological degeneracy. 

PACS numbers: 



The search for spin liquid ground state is a central issue 
in the study of the strongly correlated electron systems. 
The spin liquid state represents a novel state of matter 
beyond the Landau-Ginzburg description and supports 
new kinds of order and excitation. It is generally be- 
lieved that the study of spin liquid state will not only 
deepen our understanding on the organizing principles of 
condensed matter systems, but also result in new appli- 
cations that is not possible from conventional materials. 
In particular, the topological order and the related frac- 
tionalized excitations are proposed to implement the key 
steps of topological quantum computation. 

It is generally believed that the geometrically frus- 
trated quantum antiferromagnetic systems are ideal 
places to find spin liquid ground state. For this reason, 
the study of Heisenberg model on Kagome lattice and tri- 
angular lattice have received considerable interests. More 
recently, the spin liquid states are also proposed to ap- 
pear in Hubbard models on frustrated lattices as a result 
of the multi-spin exchange processes. The bipartite lat- 
tice, on which the antiferromagetic exchange interaction 
is not frustrated, is generally not believed to be favor- 
able for the formation of the spin liquid ground state. 
For example, on the square lattice, the system develops 
magnetic long range order as soon as one turns on the 
electron correlation. With these understandings in mind, 
it is quite unexpected that a spin liquid ground state can 
exist in the Hubbard model on the bipartite honeycomb 
lattice. 

Accordingto the numerical results reported recently by 
Meng et. al[l|, a state with full gaps to both charge and 
spin excitations and full symmetry of the Hamiltonian 
emerges in a small interaction range of 3.5 < U/t < 4.2 
for the Hubbard model on honeycomb lattice. This state 
intervenes between a semimetal phase with a Dirac-type 
spectrum for U/t < 3.5 and an antiferromagnetic ordered 
phase with spin wave excitation for U/t > 4.2, resulting 
in a counterintuitive non-monotonic evolution of the spin 
excitation spectrum with U/t. What makes this state 
even more special is that although it has a full gap to 
spin excitation, it has no accompanying topological de- 
generacy. 

In a loose sense, a spin liquid state can be defined as 



a quantum disordered insulting state of a many electron 
system that respect all the symmetries of the Hamilto- 
nian. However, to exclude the trivial case of a band in- 
sulator, in which the spin degree of freedoms cancel out 
in each unit cell, one should add the further requirement 
that each unit cell of the lattice contains an odd number 
of electrons. A system with an odd number of electrons 
in each unit cell will inevitably posses gapless excitations 
at the mean field level if there is no symmetry break- 
ing mechanism to enlarge the unit cell|3|. Built on the 
higher dimensional generalization of the celebrated Lieb- 
Schultz-Mattis theorem in one dimension, it is now gen- 
erally believed that a spin liquid state should either be 
gapless, or, while possessing a bulk gap to spin excitation, 
show topological degeneracy The topological degen- 
eracy denotes the degeneracy between ground states that 
can not be distinguished locally but are globally distinct. 
The existence of topological degeneracy is one of the most 
important manifestation of topological order, which is be- 
lieved to be a prerequisite for the emergence of fraction- 
alized excitations in the spin liquid background 4j. For 
these reasons, it is not surprising that the negative result 
on the topological degeneracy reported by Meng et.al. in 
the spin liquid state has aroused much interests and also 
confusions in the community , 5,-8] • 

In a sense, the absence of topological degeneracy in 
the Hubbard model on honeycomb lattice should not be 
that surprising since the honeycomb lattice is a complex 
lattice with two inequivalent sites in each unit cell. A 
state with one electron per site on average actually cor- 
responds to integer filling rather than half filling. Thus 
the above mentioned hypothesis on the topological de- 
generacy should not apply here. However, it is not clear 
to what extent is the oddness of electron number in each 
unit cell essential for the existence of the topological de- 
generacy in a fully gaped system. As the topological 
degeneracy is the key to the exoticness of a spin liquid[3| 
and as some doubts have been raised on the negative 
results of topological degeneracy^, II] in the numerical 
work, it is valuable to understand if Hubbard model on 
honeycomb can really support topological degeneracy. 
Even leaving apart the topological degeneracy, it is still a 
challenging task to compose a picture for the spin liquid 



state reported in the numerical works and to understand 
how a full gap can be generated in the spin excitation 
spectrum without breaking any symmetry [5|-l8|. 

In this paper, we provide an analytical argument for 
the absence of topological degeneracy in the Hubbard 
model on honeycomb lattice. This argument is based 
on the unique sign structure of the ground state of this 
bipartite system, which is shown to be insensitive to the 
trapped gauge flux when the system is defined on a torus. 
To illustrate the effect of the ground state sign structure 
on the topological degeneracy, we also present the results 
of a variational study on some Mott insulating states on 
the honeycomb lattice. 

The model studied in this paper reads 

11 = -t {c\ aCj^a + h.c.) + C/ ^ rii^n,^, (1) 

<i,j>,a i 

in which the first sum is over nearest neighboring sites 
on the honeycomb lattice. The honeycomb lattice is 
a complex lattice with two inequivalent sites in each 
unit cell(see Fig.l) and can be divided into two sub- 
lattices(even and odd). The hopping terms are nonzero 
only between sites on different sublattices. With bi and 
62 denoting the primitive lattice vectors in reciprocal 
space(6i ■ aj = 2n6ij), the momentum in the Brillouin 
zone can be parameterized as fc = fci6i + k2b2- The ki- 
netic part of the Hamiltonian can then be diagonalized 
as follows 

Ho = 5]i?M7L,.7fc<x,s, (2) 

ka.s 

in which Ek.s = s\l + e^^^ + e*'^^! and s = ±1. ^kG,s is 
corresponding eigen-operator. The bare dispersion has 
two Dirac points at (fci,fc2) = ±(27r/3, — 27r/3) in the 
first Brillouin zone. When the electron density is such 
that there is one electron per site on average, the Fermi 
surface of the free system shrinks into these two Dirac 
points. Thus, although the honeycomb lattice is bipartite 
and the antifcrromagnetic interaction is not frustrated, a 
finite strength of local correlation is needed to induced 
the instability of spin density wave ordering. 

It should be noted that the gaplessness of the bare dis- 
persion at (^1,^2) — ±(27r/3, — 27r/3) is protected by the 
bipartite nature and the three-fold rotational symmetry 
of the honeycomb lattice. Supposing that the hopping 
integral is nonzero only between sites on different sub- 
lattices, then by the three-fold rotational symmetry, a 
hopping term between sites separated by an arbitrary 
vector Ri = (m + \)ai + [n + \)o,2 is always accom- 
panied by two other hoppings of the same amplitudes 
at separations R2 = — (rn + n + |)ai + (m -I- ^02) and 
R3 = {n + ^)ai — {m + n + |)a2. It is then easy to 
check that the expression J2i=i 2 3 e*'^'^' is identically 
zero at (^1,^2) = ±(27r/3, — 27r/3). The same argu- 
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FIG. 1: The honeycomb lattice and its unit cell. 



ments can also be used to show that BCS pairing be- 
tween sites on the different sublattices and respecting the 
three-fold rotational symmetry can neither open a gap at 
(fci,fc2) = ±(2^3,-2^3). 

With this in mind, it is then quite unusual that a full 
gap in both spin and charge sector can be opened with- 
out any symmetry breaking. Could it be possible that 
the system spontaneously violates the bipartite nature of 
the honeycomb lattice by generating hopping or pairing 
terms between sites on the same sublattice[3, Q? The 
analytical argument present below suggests that this is 
very unlikely to be the case. The same argument also 
provides an understanding on the absence of the topo- 
logical degeneracy in the Hubbard model on honeycomb, 
or more generally, bipartite lattice. 

Our argument is based on the now well known Lieb's 
theorem on the sign structure of the ground state on bi- 
partite lattice for the Hubbard model. On a bipartite 
lattice such as the honeycomb lattice studied in this pa- 
per, the Hubbard model is particle-hole symmetric. To 
be more specific, through the following unitary transfor- 
mation 




in which rji = —1 for i in the odd sublattice and is oth- 
erwise 1, the Hubbard model takes the form 

il=-t i4,aC3,'^ + h-C.)-U"Yn,^n,^, (4) 

up to a chemical potential term. Under such a transfor- 
mation, a repulsive Hubbard model is mapped into an 
attractive Hubbard model. At the same time, the spin 
and charge degree of freedoms interchange their roles in 
the particle-hole transformed system. For example, a re- 
pulsive system with one electron per site is mapped into 
an attractive system with zero magnetization. 

In a celebrated paper appeared in 1989 by Lieb[l(J|, it 
is proved that the attractive Hubbard model Eq.(4) on 
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any lattice has a unique ground state with a weh defined 
sign structure. More specifically, if we use the Fock basis 
l^'t) = rifccL^^^lO) and = Uk4^u\^) ^ expand 

the ground state as = ^ 4>a,p\c(,f) <8i |/3,i), then 
the coefficient matrix </)(whose matrix elements are am- 
plitudes of the ground state in the Fock basis) can be 
shown to be positive definite, although in general its ma- 
trix elements can be negative as a result of the Fermion 
sign. This can be understood as a result of the local 
attractive interaction which encourages the down spin 
electrons to follow the trail of the up spin electrons. In 
particular, the ground state amplitude is positive definite 
when electrons of both spins share the same set of lat- 
tice sites. In term of the repulsive Hubbard model before 
the particle-hole transformation, this corresponds to a 
singly occupied configuration. Keep in mind the alterna- 
tive sign r]i in the particle- hole transformation Eq.(3), it 
can be shown easily that in the repulsive Hubbard model, 
the ground state amplitudes in the singly occupied sub- 
space satisfy the so called Marshall sign rule. The Mar- 
shall sign rule claims that the ground state amplitudes 
in the singly occupied subspace are real and their signs 
are given by (— up to a global phase factor. Here 
j\feven number of down spin electrons in the even 

sublattice. 

The existence of such a sign structure in the ground 
state, which is a result of bipartite nature of the hon- 
eycomb lattice, makes it very unlikely that the reported 
spin gap isgenerated by intra-sublattice hopping or pair- 
ing terms [J Hj. At the same time, the sign structure of 
the ground state provides a way to understand the ab- 
sence of the topological degeneracy in this system. To ap- 
preciate this point, it is crucial to realize that Lieb's proof 
is general enough to address not only the sign structure 
of the ground state, but also it response to the trapped 
gauge flux when the system is defined on a torus. Lieb's 
proof applies whenever the lattice is bipartite, while up 
to a gauge transformation trapping a gauge flux in the 
holes of the torus amounts simply to change the phase of 
the hopping integrals on the bonds across the boundary 
of the system. Thus, the ground state after the gauge 
flux trapping has exactly the sign structure as that be- 
fore the gauge flux trapping up to an inessential global 
phase. It is just this boundary condition independence 
of the sign structure of the ground state that enables us 
to arrive at our conclusion. 

To be more specific, the discussion below will be re- 
stricted to the topological degeneracy of the Z2 type, 
which is the most commonly envisaged in a gapped quan- 
tum spin liquid Q. A scheme to detect topological degen- 
eracy in the resonating valance bond(RVB) state gener- 
ated from Gutzwiller projection of BCS mean field state 
is proposed by Ivanov and Senthil in 2002 J^]. They pro- 
posed to check the orthogonality between RVB states 
with different number of trapped visons(the Z2 topo- 
logical excitation) in the holes of a multiply connected 



manifold. For a topological ordered system, states with 
different number of trapped visons in the holes are or- 
thogonal to each other and form topological degeneracy. 
On the other hand, for a topological trivial state, the 
visons can tunnel through the bulk of the system and es- 
cape from the holes. The topological degeneracy is lifted 
by such tunneling events. Built on the conceptual link 
in the RVB picture between a quantum spin liquid and 
a quantum phase disordered superconductor, Ivanov and 
Senthil also proposed to construct the vison excitation 
from quantized superconducting vortex in the mean field 
state. It is then found that as a result of the Z2 char- 
acter of the quantized superconducting vortex, trapping 
a vison in a hole of multiply connected manifold for the 
RVB state amounts to changing the boundary condition 
around the hole from periodic to anti-periodic or vice 
versa. 

In a later generalization by the present author and 
Yang [13], the role of the superconducting vortex in defin- 
ing the vison is replaced by that of a Z2 gauge flux. Thus, 
to detect topological degeneracy in the RVB state, one 
should check the sensitivity of the ground state to a quan- 
tized gauge flux in the holes on a multiply connected 
manifold. Such a deflnition closely resembles the way 
that Kohn chose to define a insulator ll| . According to 
this definition, the key difference between an insulator 
showing topological degeneracy(order) and a trivial in- 
sulator is that the former can respond nonlocally to a 
trapped Z2 gauge flux while the latter has no such non- 
local response. 

In 13| , it is also shown that RVB states generated from 
bipartite mean fleld Hamiltonian through Gutzwiller pro- 
jection all satisfy the Marshall sign rule, no matter what 
is the gauge structure of the corresponding effective the- 
ories at the Gaussian level. As a result of such sign struc- 
ture, it is proved generally that all RVB state generated 
from a bipartite mean field Hamiltonian can not sup- 
port topological order. The same argument can also be 
adopted here after some modifications. The modification 
is necessary for two reasons. First, here we are discussing 
ground state of a given Hamiltonian rather than a gen- 
eral RVB state with unspecified Hamiltonian. Second, 
the model we discuss here also has charge excitation, the 
existence of which will inevitably cause Fermion sign and 
the ground state wave function has a definite sign only 
in the subspace of singly occupied state. 

To deal with the first problem, we simply introduce the 
Z2 gauge flux directly in the Hamiltonian rather than in 
the effective theory as was done in the previous study. 
From the above argument on the topological degeneracy, 
this seems to be intuitively even more appealing as it 
allow us to probe directly the response of the ground 
state to the trapped Z2 gauge flux. 

Now we check the orthogonality of the ground states 
with and without a trapped Z2 gauge flux in the holes of 
a multiply connected manifold. Let (/) and be the wave 
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function matrixes of the ground states \^) and j^" ) before 
and after the flux trapping, then as the sign structure of 
the ground state is insensitive to the trapped gauge flux, 
both matrixes should be positive definite up to a global 
phase factor. The overlap of the two states is given by 
(^'1^' ) = Ticfxj) , in which we have used the fact that 
both (p and (f> can be assumed to be Hermitian without 
loss of generality. In the diagonal representation of (j), the 
overlap can be written as ) = J2i'^i4'iiJ i^i which 

Wi are the positive definite eigenvalues of the matrix (jj 
and (j)^ ^ are the diagonal matrix elements of matrix (p 

is this representation. Since is also positive definite, 
0j j > 0. Thus the overlap is the sum of positive terms 
and is nonzero for finite system. 

To prove the absence of the topological degeneracy, 
one should still go to the thermodynamic limit to see 
if the overlap converges to a finite value. At present a 
proof of this is beyond our reach. With this in mind, 
the results presented above can be either interpreted as 
a suggestive evidence for the absence of the topological 
degeneracy, or, on the other way around, an excuse for 
the failure to detect the anticipated topological degener- 
acy on finite system. However, forbiddenness of topolog- 
ical degeneracy induced by the Marshall-sign structure 
of ground state has been illustrated before at the vari- 
ational level for various RVB states and the same thing 
can also be done for the honeycomb system. 

For this purpose, we investigate two typical Mott in- 
sulating states on the honeycomb lattice with built-in 
Marshall sign structure. Both of these states are gener- 
ated by Gutzwiller projection of BCS mean field states. 
The first is the Gutzwiller projection of the free Fermion 
state on the honeycomb lattice with nearest neighbor- 
ing hopping terms and two Dirac points at (fci,fc2) — 
±(27r/3, — 27r/3). The second is the Gutzwiller projec- 
tion of the BCS pairing state with d + idl pairing pattern 
between neighboring sites lj| . The Z2 gauge flux can be 
imposed in the boundary condition of the BCS mean field 
Hamiltonian[ll[l3. 

As a result of a general theorem proved in [l^, both 
states satisfy the Marshall sign rule, even if the unpro- 
jected wave function is not real(as for the d + id' state). 
To see if topological degeneracy can survive in these 
states, we calculate the overlap of the states defined on a 
torus with different number of Z2 gauge flux threaded in 
both holes of the torus. The overlap between the wave 
functions can be easily calculated by with the variational 
Monte Carlo method. The results for the projected free 
Fermion state is presented in Fig. 2, in which we have 
shown the overlap between the state with no flux in both 
holes and the state with Z2 flux in one hole and no flux 
in the other hole of the torus. The calculation is done on 
lattice with Lx Lx2 sites and L is so chosen to avoid the 
Dirac nodes in the momentum space. The results rapidly 
converges to a finite value when L > 10 . 
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FIG. 2: The overlap between states with different number of 
Z2 gauge flux in both holes of the torus for the Gutzwiller 
projected free Fermion state on honeycomb lattice. I^*) de- 
notes the state with no gauge flux in both holes of the torus 
and Ivp') denotes the state with Z2 gauge flux in one hole but 
no flux in the other hole of the torus. 



An even more striking example of how the sign struc- 
ture of the ground state can cause the absence of the 
topological degeneracy is provided by the projected d+id' 
wave pairing state on the honeycomb lattice. Here, 
the unprojected mean field state break the time rever- 
sal symmetry and has a complex wave function. In 
terms of the effective field theories joj, the RVB state 
generated from such a mean field ansatz would posses 
a Z2 gauge structure and thus support topological or- 
der. However, after the Gutzwiller projection , the wave 
function becomes real up to a global phase and satisfy 
the Marshall sign rule. In Fig. 3 we present the over- 
lap of the state with Z2 gauge flux in the hole sur- 
round by the x-circumference but no flux in the hole sur- 
rounded by the y-circumference(denoted as l^'io)) and 
the state with Z2 gauge fiux in the hole surround by the 
y-circumference but no flux in the hole surrounded by 
the x-circumference (denoted as j^'oi)). The overlap is 
also seen to converge to a finite value in the thermody- 
namic limit. 

Although both variational states presented above are 
not meant to describe the spin liquid state found in the 
numerical simulation of [l[, they do illustrate well how 
the sign structure of the ground state can cause the non- 
orthogonality of the states with and without trapped 
gauge flux in the thermodynamic limit and thus the for- 
biddenness of the topological degeneracy of the system. 
A full description of the spin liquid state found in [l[ is 
beyond the scope of this paper and requires obviously 
further refinement of the variational wave function. For 
example, the virtual charge fluctuation should obviously 
be taken into account. However, we think these modifi- 
cations will not change our conclusion qualitatively. 
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FIG. 3: The overlap between states with different number of 
Z2 gauge flux in both holes of the torus for the Gutzwiller 
projected d + id! pairing state on honeycomb lattice. |^'io) 
denotes the state Z2 gauge flux in the hole surround by the 
x-circumference but no flux in the hole surrounded by the y- 
circumference and I'&oi) denotes the state with Z2 gauge flux 
in the hole surround by the y-circumference but no flux in the 
hole surrounded by the x-circumference. 

In summary, we have presented both analytical argu- 
ments and variational examples that suggest the topolog- 
ical degeneracy can not survive in the Hubbard model on 
honeycomb lattice as a result of the unique sign structure 
of its ground state. The same sign structure is believed 
to be important also for other aspects of the ground state 
and should be taken into account in the variational study. 
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